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ON STABLE COMPACTIFICATION WITH CASIMIR-LIKE POTENTIAL 
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Multidimensional cosmological models with a higher dimensional space-time manifold M — R x Mo x Y17-1 Mi 
(n > 1) are investigated under dimensional reduction to Do -dimensional effective models. In the Einstein conformal 
frame, the effective potential for the internal scale factors is obtained. The stable compactification of the internal 
spaces is achieved due to the Casimir effect. In the case of more than one internal space a Casimir-like ansatz for 
the energy density of the massless scalar field fluctuations is proposed. Stable configurations with respect to the 
internal scale factor excitations are found in the cases of one and two internal spaces. 
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1. Introduction 

Our observable universe at the present time at large 
scales is well described by the Friedmann model with 
4-dimensional Friedmann-Robertson- Walker metric. 
However, it is possible that space-time at short (Planck) 
distances might have a dimensionality of more than 
four and possess a rather complex topology M. Ac- 
counting this possibility it is natural to generalize the 
Friedmann model to a multidimensional cosmological 
model (MCM) with the topology B, M 



M = IxIoxIiX...x M n , 



(1.1) 



where M, (i = 0, . . . , n) are spaces of constant cur- 
vature (more generally, Einstein spaces). Mq usually 
denotes a do -dimensional external space and Mi (i = 
1, . . . , n) are di -dimensional internal spaces. These in- 
ternal spaces have to be compact, what can be achieved 
by appropriate periodicity conditions for the coordi- 
nates (§-]§]■ As a result, internal spaces may have 
a nontrivial topology, being compact (i.e. closed and 
bounded) for any sign of the spatial curvature. 

To make the internal dimensions unobservable at 
the present time the internal spaces have to be reduced 
to scales near the Planck length Lpi ~ 10~ 33 cm, i.e. 
scale factors of the internal spaces Oj should be of or- 
der of Lpi . Obviously, such compactifications have 
to be stable. It means an effective potential of the 
model obtained under dimensional reduction to a 4- 
dimensional effective theory should have minima at 
aj ~ Lpi (i = 1, . . . ,n). A lot of papers were de- 
voted to the problem of stable compactification of ex- 
tra dimensions (an extended list of references may be 
found in the paper ||). A number of effective poten- 
tials which ensure stability was obtained there. Among 
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them the Casimir potential is one of the most impor- 
tant. The Casimir effect is connected with the vac- 
uum polarization of quantized fields due to non-trivial 
topology of the background space or the presence of 
boundaries in the space. As a result one obtains a 
nonvanishing energy density of the quantized fields in 
the vacuum state. In our case this phenomenon should 
take place due to internal spaces compactness. 

The Casimir energy density for a massless scalar 
field in a model with one finite scale factor has a form 

ifuj-iEa 



(1.2) 



where D is the total space-time dimension and C is 
a constant which depends strongly on the topology of 
the model. Equation (1.2) holds for models with one 
factor-space (e.g. Mq with the scale factor ao = a) 
@, ||, @, 0], || , or in the case of one internal 
space and an external scale factor that is much greater 
than the internal one: ao di = a ]l2[-Jl4[. 

In the case of several internal compact spaces the 
Casimir effect occurs if at least one of them is finite. 
It is natural to consider all internal spaces on equal 
footing supposing that their scale factors are freezed 
near Planck length. In the papers |2^, [24j] an approxi- 
mation of the Casimir energy density was proposed in 
the form 



f 



n 



l — 1 l—l 



am 

„4 



(1.3) 



for the universe with the topology M = 
n™=i 3 d where Vi ~ af is the volume of the d - di- 
mensional sphere. It can be easily seen from this equa- 
tion that p — ► if any of a$ — > oo . This means that 
approximation (|0|) is not applicable in this limit be- 
cause, first, for any finite a^ there should exist a nonva- 
nishing Casimir energy density, and, seco nd, it should 
provide a correct transition to Eq. (1.2), what, obvi- 
ously, is not the case. 
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The calculation of the Casimir energy density even 
in the case of one scale factor is a complicated problem. 
For several scale factors this procedure becomes ex- 
tremely difficult, especially analytically. In the present 
paper we suggest ad hoc a Casimir-like ansatz for the 
energy density of the massless scalar field fluctuations, 
which yields a better approximation than flL3| ). In this 
ansatz all internal scale factors are included on equal 
footing, it has the correct Casimir dimension: [cml D 



and gives a correct transition to the formula (1.2). The 
proposed ansatz gives an energy density that does not 
equal to zero if at least one of scale factors is finite. 

As it was mentioned above the problem of stable 
compactification is one of the most important in MCM. 
In our paper we investigate this problem for the pro- 
posed Casimir-like potential in the case of one and two 
internal scale factors. It is shown that stable configu- 
rations exist for both of these cases. To our knowledge 
the investigation of the stable compactification prob- 
lem due to Casimir potential for models with more 
than one non-identical (with respect to the scale fac- 
tors) spaces was not performed up to now. 

The paper is organized as follows. In Section 2, the 
general description of the considered model is given. 
In Section 3, the effective potential is obtained under 
dimensional reduction to a Dq -dimensional (usually 
Dq = 4) effective theory in the Einstein frame. The 
problem of stable compactification is investigated in 
Section 4 for one internal space and in the Section 5 
for two internal spaces. Stable configurations are found 
for both of these cases. Some arguments in favour of 
our Casimir-like ansatz are listed in the Appendix A. 
In Appendix B, some general features of a generaliza- 
tion of the Abel-Plana summation formula to higher 
dimensional complex spaces are described. A simple 
example for a potential with nondegenerate minimum 
at a point with coinciding scale factors is given in Ap- 
pendix C. 

2. The model 

Let us consider a cosmological model with the metric 



ds 1 



-e^dr 2 



(2.1) 



i=0 



which is defined on the manifold (1.1), where the man- 
ifolds Mi with the metrics g^ are Einstein spaces of 
dimension di , i.e. 



Rr, 

and 
R 



(■') 



aV i) 



0, . . . ,n 



— A'cL = Rj 



(2.2) 



(2.3) 



In the case of constant curvature spaces X 1 = ki(di — 1), 
ki — ±1,0. The non-zero components of the Ricci 



tensor for the metric fl2.1|) read 

n r 

i?oo = -X)di /3 l -7/3+(/3 1 )' 



R,, 



i=0 



(2.4) 



g£\ nt [A l + exp(2/3 l -27) X 
x UT+p (f^d^-jj) , (2-5) 
so that the corresponding scalar curvature is given as 

n n 

R = ^i? 2 exp(-2/3 l ) +exp(-2 7 )^d 4 x 



i=0 



i=0 



3=0 



(2.6) 



The overdot in formulae ( p.4| ) - (2J3) denotes differen- 
tiation with respect to time t . 

The action of the model we choose in the following 
form 

S =^J d D x^\g~\{R[g}~2A} + S c + S YGH} (2.7) 

M 

where D = 1 + J^ILo d% * s ^ ne total dimension of the 
space-time, A is a D -dimensional cosmological con- 
stant, k 2 — D -dimensional gravitational constant and 
Sygh is the standard York-Gibbons-Hawking bound- 
ary term. The Casimir effect is taken into account via 
the additional term S r of the form |25l 



S c = - J d D x-y/\g\p(x), 

M 



(2.8) 



where p{x) is the Casimir energy density (A.ll) . For 
the model under consideration the Casimir energy is 
a function of the scale factors (see ( A. 28]) - (A. 30)). 
After dimensional reduction the action (2.7) reads 



5' 



drL 



(2.9) 



with Lagrangian L given as 



n n 



-2/3- 



i,j=0 

K 2 

e 7+70 A e 7 Fc _ 

M 



i=0 



(2.10) 



Here 70 = Yh=o an d f£ = Il"=o Mi where fa is 
defined by the equation ( A.12[ ). The components of 
the minisuperspace metric read (3) 



(^ij d^S^j d^ dj , % , j , . 



(2-11) 



It is easy to show tha t the Euler-Lagrange equa- 
tions for the Lagrangian ( 2.1 0| ) are e quivalent to the 
Einstein equations for the metric (|2.l|) and the Casimir 
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energy-momentum tensor ( A. 14 ). This follows from 
the fact that for the non-trivial components of the ten- 
sor Emn we get 



E, 



(HI 



dL 



exp (7-70), 



1 w f d_dL 



dL 

x exp (2(3 l - 7 - 70), 
i = 0, ... ,n , 

where the tensor Emn is defined as 
1 



(2.12) 



(2.13) 



E 



MN 



Rmn — -z9mnR ~ k T j 



MN 



A 5 



MN- 



(2.14) 



For t he pr oof on e makes explicit use of equations 
(|A 11|) and flA.13|) . 



3. The effective potential 

Let us slightly generalize this model to the inhomo- 
geneous case supposing that the scale factors (3 l — 
f3 l (x) (i = 0,...,n) are functions of the coordinates 
x , where x are defined on the Dq = (1 + do) - dimen- 
sional manifold M = Ix M with the metric 

= gf) dx » g, dx» = -e^dr 2 + e 2 ^y°). (3.1) 

Mq denotes the external space-time. The dimensional 
reduction of the action (^) yields 



S = ^ I d D °x^\\[e^ {R 



M 



i=l 



-,(0) 



i=l 

„2 , 



2A - 2k 2 p} 



(3.2) 



where Kq = K 2 //i is the Do -dimensional gravitational 
constant, p = 11"=!^ anc ^ ^ij = l,...,n) is 



the midisuperspace metric with the components (2.11). 
Here the internal scale factors play the role of scalar 
fields (dilatons). The action ([T^) is written in the 
Brans-Dicke frame. Conformal transformation to the 
Einstein frame 



5(0) 

it fiv 



yields 

S - 



g(0) 



(3.3) 



2kq 



d D °x^\g(°)\{R 



M 



Gijg^d^ d v fi - 2U eff } , (3.4) 



where the tensor components of the midisuperspace 
1 



metric Gy are 



Gij d^Sij 



: didj , 



Dq _2~>~>> hJ 
and the effective potential U e ff reads 



(3.5) 



U, 



<ff 



1 ™ 

-Je^ 



-2/3' 



+ A + n 2 p] . (3.6) 

In what follows we consider the case where the Casimir 
energy density depends on the internal scale factors: 
p = p (/3 1 , . . . , /3") . This means that either the ex- 
ternal space Mq is non-compact or the external scale 
factor is much greater than the internal ones: ao 3> 
01, . . . , a n ■ 

In the case of one internal space n = 1 the action 
and the effective potential are respectively 

* - 4i / d D ° x J\g~w\{R 



2n r . 



M 



g^^d^d^-2U eff } 



(3.7) 



and 

U e ff = e 2tP [(D-2)'( 1 Do-2)] 



„ Dp-2 



where we redefined the dilaton field as 



l dl (D 



D 



(3.8) 



(3.9) 



The minima of the effective potentials ( |3.6| ) and ( p.8| ) 
define stable compactification positions = exp (3 l c , 
c = l,...,m and small excitations ip l of the inter- 
nal scale factors near these positions have a form of 
minimally coupled massive scalar fields in the external 
space-time Q: 

1 



S 



2kq 



d Do x^\g(°)\{R 



(0) 



2A 



(c)e// 



} 



Mo 



n _ 

d D °x^\g(°)\{-g^^- 



M u 



(3.10) 



where A-( c ) e ff = U e ff(f3 c ) is an effective cosmological 
constant and m( c )^ are masses of gravitational excitons 
corresponding to the c-th minimum. From a physi- 
cal point of view it is clear that the effective potential 
should satisfy following conditions: 



(<) 
(ii) 

(in) A (c ) e// 



«(c)i 
™(c)i 



< 



Lpi , 
. 
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The first condition expresses the fact that the inter- 
nal spaces should be unobservable at the present time 
and stable against quantum gravitational fluctuations. 
This condition ensures the applicability of the classi- 
cal gravitational equations near positions of minima of 
the effective potential. The second condition means 
that the curvature of the effective potential should be 
less than Planckian one. Of course, gravitational exci- 
tons can be excited at the present time if vrn <C Mpi . 
The third condition reflects the fact that the cosmo- 
logical constant at the present time is very small: A < 
lCT 56 cni- 2 w lCT 121 Api, where A P i = L p 2 . Thus, 
for simplicity, we can demand A e // = U e ff((3 c ) = 0. 
(We used the abbreviation A e // = Afcjeff ■) Strictly 
speaking, in the multi- minimum case (c > 1) we can 
demand d( c )j ~ Lpi and A/ C \ e ff = only for one of 
the minima to which corresponds the present universe 
state. For all other minima it may be ei( G )j 3> Lpi and 



A, 



c)eff 



>0. 



4. One internal space 

Here we consider the case of one internal space n = 1 
or, stricly speaking, the case where all internal spaces 
have one common scale factor: a\ = a® = . . . = a n = 
a = cxp fj . The latter case corresponds to the split- 
ting of the internal space into a product of Einstein 
spaces JM: Mi — > 11™= i % which leads in the equa- 
tions ( |3.7|) - ( fTs| ) to the substitutions: d% — ► d = 
2™=i 4 i Ri —* R = Yl?=i Ri ■ For effective one- 
scale-factor models the Casimir energy density reads 
according to ( A. 26 ): p = Cexp(— D/3), where C is a 
constant that strongly depends on the topology of the 
model. For example, for fluctuations of massless scalar 
fields the constant C was calculated to take the values: 
C = -8.047- 10- 6 if M = Rx S 3 , M x = S 1 (with e fJ ° 
as scale factor of S 3 and e 13 " > ) Jlj 



if M 
if M 



X K 
X S 3 



Mi = S 1 JT| and C 
Mi = S 3 (with e f3 ° > e 



C = -1.097 
: 3.834- 10- 6 





For an effective potential with p 



the zero-extremum-conditions 



d/3 



Cexp(-L> / 3) 
= and 



A e // = lead to a fine tuning of the parameters 
of the model 



Hie 



-20c 



2L> 
£>-2 



A, i? ie ( D ~ 2 )^ = k 2 CD, (4.1) 



which implies signi?i = sign A = signC. We note 
that a similar fine tuning was obtained by different 
methods in papers (for one internal space) and |27[ ] 
(for n identical internal spaces). To get a minimum, 
the second derivative of the effective potential at the 
extremum should be positive: 



d 2 U, 



eft 



dp 2 



= (D- 2)i?i (e- 2 ^) D °- 2 > 0. 



(4.2) 



Thus, a stable compactification takes place if the in- 
ternal space has a positive scalar curvature R\ > (or 
for a split space Mi the sum of the curvatures of the 
constituent spaces M* should be positive). The mass 
squared of the gravitational excitons is given by the 
expression 



d 2 u t 



eff 



dip 2 



Do 



di 



-Rx (e- 2 ^) D °- a . (4.3) 



As example, let us consider a manifold M with 
topology M = R x S 3 x S 3 where a > a±. Then 
according to |Q the constant C is given as C = 
3.834 • 10" 6 > 0, so that with C,Rx = di(d x - 1) = 
6 > the effective potential has a minimum provided 
A > 0. Normalizing Kq to unity, we get n 2 = p, where 
p = 27r( d+1 )/ 2 /r (i(d+ 1)) is the volume of the d- 
dimensional sphere. For the model with the topology 
M = Rx S 3 x S 3 we obtain a c m 1.5 • 10 _1 L P; and 
m w 2.12 • 10 2 Mp; . Hence, for this particular example 
the conditions (i) and (ii) are not satisfied . It would 
be interesting to investigate models with more complex 
topology. Of course, the calculation of the Casimir ef- 
fect in this case becomes a very complicated problem. 
Therefore in the Appendix A we propose ad hoc the 
Casimir-like expressions ( |A.28|) , ( |A.29| ) for the energy 
density in the case of more than one internal spaces 
with non- identical scale factors. In the next section 
we investigate the stable compactification problem for 
this potential in the case of two internal spaces. For 
this case the Casimir-like potential considerably sim- 
plifies, but even in this case the stability analysis is 
still complicated. To our knowledge the investigation 
of the stable compactification problem for models with 
more than one non-identical (with respect to the scale 
factors) spaces was not performed up to now. 

5. Two internal spaces 

After these brief considerations on Casimir potentials 
for one-scale-factor models we turn now to some meth- 
ods applicable for an analysis of two-scale-factor mod- 
els with Casimir-like potentials. Some arguments in 
favour of such potentials are given in Appendix A. We 
propose to use these potentials of the general form 

fe o ,...,kn=0 £o=0£i=0 



(e^f...^")^ 1 



V 4 (feo) 

?0...«n-l / j gfc \fo+€l + -+f„_l + l ' 



£»-i=o 



(5.1) 



in order to achieve a first crude insight into a possi- 
ble stabilization mechanism of internal space configu- 
rations due to exact Casimir potentials depending on 



On stable compactification with Casimir-like potential 



5 



n scale factors. In ( |5.l| ) A^ ^ _ i are dimensionless 
constants which depend on the topology of the model 
and e lk ... m = (-l) P £;fc...m ■ Here e ik ... m is the to- 
tally antisymmetric symbol (eoi...m = +1), P is the 
number of the permutations of the 01 ... n resulting in 
ik . . . m . \k1k2 ■ ■ ■ k n \ means summation is taken over 
k\ < k 2 < ■ ■ ■ < k n . 

From investigations performed in the last decades 
(see e.g. [jl6|, [l7|, Rcfs. therein and Appendix A be- 
low) we know that exact Casimir potentials can be 
expressed in terms of Epstein zeta function series with 
scale factors as parameters. Unfortunately, the exist- 
ing representations of these zeta function series are not 
well suited for a stability analysis of the effective po- 
tential U e f f as function over the total target space (3 £ 

. The problems can be circumvented partially by 
the use of asymptotic expansions of the zeta function 
series in terms of elementary functions for special sub- 
domains Q, a of the target space il a C M.^ . According 
to [jl6|, |l7j potential ([[j]) gives a crude approximation 
of exact Casimir potentials in subdomains fl a ■ In con- 
trast with other approximative potentials proposed in 
literature |23|, potential (5.1) shows a physically 
correct behavior under decompactification of factor- 
space components. The question, in as far ( |5.lD can be 
used in regions \ fl a , needs an additional investiga- 
tion. The philosophy of the proposed method consists 



in a consideration of potentials (5.1) on the whole tar- 
get space > an( i testing of scale factors and param- 
eters of possible minima of the corresponding effective 
potential on their compatibility with asymptotic ap- 
proximations of exact Casimir potentials in Q a . As a 
beginning, we describe in the following only some tech- 
niques, without explicit calculation and estimation of 
exciton masses. 

Before we start our analysis of two-scale-factor 
models with Casimir-like potentials 



■>(D +i)P 



3=0 



(2) e^ 1 

3 e (A>+j)/3 2 



(5.2) 



let us introduce the following convenient (temporary) 
notations: x := a\ = e@ , y := a% = eP , := 

the effective potential (|3.6|) reads 



k 2 A?\ In terms of these notations 



U, 



eff 



(x dl y d2 ) 



Ri - 



<i 2 



+ x - dl y- d2 \ Y^P l y % x 

\i=0 



2 



3=0 



3 x 3 y 



<D +j) 



(5.3) 



For physically relevant configurations with scale factors 
near Planck length 



< x, y < 00 



(5.4) 



we transform extremum conditions dpiaU e ff = 
d x ,yU e ff = by factoring out of (xy)~ D — terms 
and taking combinations d x U e ff ± d y U e ff — to an 
equivalent system of two algebraic equations in x and 
V ■ 

D- 2 , , 

(Riy 2 



h+ = (*y) D - 2 

2A 



R 2 x 2 



)- 



D x 2 y 2 



2D 



Dn 



D x 



J2Piy D " +dl+ixd2 ~ l 

i=Q 



j=0 



S S V' 



h-j x D +d 2 +j 



= 



h- = (xy) 



D-2 



di 



Dn 



+ (R iy 2 - R 2 x 2 ) - 



2 v 
2A 



3=0 
0. 



D () 



Do 

di - d 2 
D -2^ 

d\- d2 _ 
D -2 ~ 



Riy 2 + R 2 x 2 



(d 1 -d 2 )x 2 y 2 



2/ 



2j 



yDo+di+i x d2—i 



ydi—j x D +d 2 +j 



(5.5) 



(Here we used the notation D := d\ + d2-) Thus, scale 
factors ai and (Z2 satisfying the extremum conditions 



are defined as common roots of polynomials (5.5). In 
the general case of arbitrary dimensions (Z?0j di, ^2) 
and arbitrary parameters {R\, R2, Pi, Si} these roots 
are complex, so that only a restricted subclass of them 



are real and fulfill condition (5.4). In the following 
we derive necessary conditions on the parameter set 



guaranteeing the existence of real roots satisfying (5.4). 
The analysis could be carried out using resultant tech- 
niques Jl8[ | on variables x, y directly. The structure of 
Ii± suggests another, more convenient method [ p"9| . 
Introducing the projective coordinate X = y/x we 
rewrite (|5.5|) as I\± — x D l2±(y, A) with 



2+ 



-a (X) + a_D_ 2 (A)y 



D-2 



h- = -60(A) + b D - 2 {X)y 
and coefficient-functions 



D-2 



y D ^ + 
y D A-- 



-0 (a) 
(b) 



(5.6) 



ao (A) 



l Do -2 



D] 



j2^ Do+di+i 

i=0 



E S 3^ 

3=0 



ac>-2(A) 



D 
Do~ 



(R 1 X 2 + R 2 



(i 
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a-1 r 
i=0 L 


^0 ( 


' ' di — d 2 




J=0 


£>0 


fdi- d 2 


\Do-2 



2i 



X 



>D-1 



A- 



(A) 



di 



Dn-2 



(RiX 2 + R 2 ) 



-2j 
'i?iA 2 



i?2 



2A 



;(di±d 2 ). 



(5.7) 



Equations ( pjq ) have common roots if the coefficient 
functions {<Zi(A), £>i(A)} are connected by a constraint. 
This constraint is given by the vanishing resultant 



R y [h+ , h 



w{\) = 0. 



(5. 



Now, the roots can be obtained in two steps. First, 
one finds the set of roots {Xi} of the polynomial w(X) . 
Physical condition (5.4) on the affine coordinates (x, y) 



implies here a corresponding condition on the projec- 
tive coordinate X — y/x 



Im(X) = 0, 



< A < 



(5.9) 



Second, one searches for each Xi solutions {yij} of 
( [5.6] ) . The complete set of physically relevant solutions 
of system ( |5.5| ) is then given in terms of pairs {xij — 
Vij/\, .'/,., I ■ 

Because of the simple y— structure of equations 
( |5.(| ) the polynomial w[X) can be derived from ( |5.6| ) 
directly, without explicit calculation of the resultant. 
Taking b (\)I 2+ -a (\)I 2 _ = 0, A_/ 2+ - A + / 2 _ = 
and assuming y > we get 



.-D-2 



Li 
L 2 



(5.10) 



where 



L X {X) 
L 2 (X) 
L 3 (X) 



(5.11) 



A_a (A)- A+6 (A) 
A^a D _ 2 {X)- A + b D _ 2 {X) 
a (A)&D_ 2 (A) - b (X)a D _ 2 (X) 

depend only on A. Excluding y from ( |5.10|) yields 
the necessary con stra int for the coefficient functions of 
equation system ( |5.6| ) 

Lt(X)Lf-*(X)-L?(X) 



w(X) 



0. 



(5.12) 



Together with condition (5.9), this polynomial of de- 
gree 

deg x [w(X)}=D 2 (5.13) 

can be used for a first test of internal space configura- 
tions on stability of their compactification. If the cor- 
responding parameters {R\, R 2 , Pi, Si} allow the exis- 
tence of positive real roots Xi , the space configuration 
is a possible candidate for a stable compactified con- 
figuration and can be further tested on the existence 



of minima of the effective potential U e f / . Otherwise it 
belongs to the class of unstable internal space configu- 
rations. 

Before we turn to the consideration of two-scale- 
factor models with factor-spaces of the same topolog- 
ical type (Mi — M 2 ) we note that for the coefficient 
functions ( ^5 .11 ), because of (5.4) and (5.1C), there 
must hold 



sign(Li 



sign(L 2 



I A; 



sign(L 3 )| 



(5.14) 



Furthermore we see from ( |5.12D that for even dimen- 
sions D = dim(Mi) + dim(M 2 ) + dim(M ) + 1 of the 
product-manifold the polynomial w(X) factors into two 
subpolynomials of degree D 2 /2 



w(X) 



x 

= 0. 



L 2 {X)L^ (\)+L?(X) 
L 2 (X)L^ {X) - if (X) 



(5.15) 



5.1. Two identical internal factor-spaces 

In the case of identical internal factor-spaces Mi and 
M 2 we have d\ = d 2 , Pi = Si, R\ = R 2 . If we as- 
sume additionally an external space-time Mq with 
dim Mp = 4 and, hence, D — 2 (d 1 + 2), then equa- 
tions (5.6) and polynomial ( |5.15 ) can be rewritten as 



-4(di + l)ao(A) + (di + l)i?i(A 2 + l)y 



D-2 



-2d x ky D = 



(a) 



(5.16) 



h- = (A 2 - 1) [-26o(A) + Riy°- 2 ] = (b) 



and 



;(A) = 4A 2 d 2 [2(A 2 - 1)] 



2(di+2) 



R dl+2 (2(di + 1)) 



di+1 ^di + l 



(-2Adi) 



°0 



w+(X) 

R dl+2 (2(di + l)) dl+1 L dl+1 - (-2Adr) rfl+1 5* 



W-(\) 

= 

with the notations 

L z := 2a -(A 2 + l)5 

a := J2 P * [X 4+dl+l + X dl - 1 ] 



(5.17) 



=o 



i+1 



b := ^P l (2 + i )A dl -^A^. 

i=0 j=0 



(5.18) 
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According to t he fac torizing polynomial ( 5.16 (b)) and 
the constraint (|5.17|) the root set splits into two types 
of subsets defined by the conditions 

type I: 

I 2+ (A, y) = 0, 7 2 _ ( A, y) = 0, w+ (A)tu_ (A) = 
type II: 

/ 2 +(A,y)=0, A 2 -1 = 0. (5.19) 
For roots of type I we have similar to ( 5.1 0| ) 
(d 1 + l)R 1 L 3 D _ 2 2b 



2Adife 



y 



Ri 



(5.20) 



whereas type II roots at A — 1 should be found from 
the polynomial I 2+ ( 5.16 (a)) directly. For this poly- 
nomial we have now simply 

/ 2+ (A = 1) := -j^y 2idl+2) - Riy 2(dl+1) + 4^> 



0. 



(5.21) 



Coming back to the general case of identical factor- 
spaces Mi, M 2 with coinciding or noncoinciding scale 
factors we note that there exists an interchange sym- 
metry between M\ and M 2 , which becomes apparent 
in the root structure of the polynomial w(X) . From 



Ueff = (xy) 



— ^-(x -y ) 



A- 



(xy)~ 



J2 Pi {y l ^~ l + x'y-^) 



i=0 



(5.22) 



we see that x and y en ter (5.22) symmetrically. When 
one extremum of ( 5.22 ) is located at {xi = a, y% = b} 
then because of the interchange symmetry x ^ y there 
exists a second extremum located at {xj = b, yj = a} . 
So we have for the corresponding projective coordi- 
nates : 



Aj = Vi/xi = b/a, Xj = y 3 /xj = a/b 
=► Ai = A7 1 . 



(5.23) 



By regrouping of terms in (5.17) it is easy to show that 
wiX- 1 ) = \- D ' 2 w(\) (5.24) 

and, hence, roots {A; ^ 0} of w(X) = exist indeed 
m pairs {Xi, A^ 1 }. But there is no relation connecting 
this root-structure with a symmetry between w + (X) 
and W-(X) in ( |5.17 ) w+(X^ 1 ) 7^ u >_ (A) . For complete- 
ness, we note that relation ( 5.23| ) is formally similar 
to dualities recently investigated in superstring theory 

a 

Before we turn to an analysis of minimum con- 
ditions for effective potentials U e ff corresponding to 



special classes of solutions of w(X) — we rewrite the 
necessary second derivatives 



Ri 



x u eff 



( ai x~ dl - 4 y- dl + 



a 2 x^- 2 y- dl - 2 ) + 
Aa 2 x- dl - 2 y- dl + 

Y^P^a^-^x- 1 - 2 ^ 



r- 

yy 





i=0 






+ a^- 2d ^ 2 y-'- 2d ^) 


Ueff 


= d lx U eff 


x^y 


Ueff 


Ri 
—as 


[x-^y-^- 1 + 




+ x~ dl ~ 1 y 






+ Aa 6 x~ dl 






di 






+ a 4 ^- 2dl - 1 y- l - 2d1 - 5 ) , 



(5.25) 



where 

a x = (d x + 2)(di + 3) 

a 2 = di(di + 1) 

a 3 = (2d 1 +i + 4)(2d 1 +i + 5) 

a± = (2di - i)(2di - i + 1) 

a 5 = dx(dx + 2) 

a 6 = d\ 

a 7 = (2d 1 +i + 4)(2di-i), (5.26) 
in the more appropriate form (notation Jl — \ dl y^ 2D+2 ) 



d 2 U, 



x u eff 



X 2 J1 



di 



~ (a 1 X 2 + a 2 )y D 2 + Aa 2 y 1 



^Pi(a 3 A 4+<ll+< + a 4 A dl - 



»=0 



dyyUeff 



Ri 



ai+a 2 X 2 )y D 2 + Aa 2 y 



,D 



d lyUejf 



^P i (a 4 A 4+dl + i + a 3 A rfl - i ) 
Ri 



i=0 



Xfi 

di 



-a 5 (X 2 + l)y D 2 + Aa 6 y 



^P 4 a 7 (A 4+dl+i + A dl - 



i=0 



(5.27) 
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Introducing the notations 

A c := 
and 



dlxUeff dlyU eff 



d ly U eff 



11' 



(c)l,2 



~ Tr(A c )± ^Tr 2 (A c )~ 4 det(A c ) 
the minimum conditions are given as 

10( C )1 > 0, W( c)2 > 0. 



(5.28) 



(5.29) 



(5.30) 



In the degenerate case of coinciding scale factors x — 
y, A = 1 there hold the following relations between 
the derivatives of effective potentials U e ff(x,y) and 
Ueff(y) = U eff {y,y) 



d y u eff 
>m" e ff 



= d x U eff \ x=y + d v U eff \ x=y 



dLu, 



d X x U eff\ x =y+ d yy U eff\ x =y 
+ 2d xyUeff\ x=y 

and minimum conditions reduce to 



dyUeff = 0, 
with 



dyyUeff 



dtyUeff > 



Adi(2di + l)y D 



(5.31) 



(5.32) 



l)(2di + 3)» 



D-2 



4(di + l)(4di + 5)^Pi 

i=0 



(5.33) 



For convenience of the additional explicit calculations 
of constraint U e ff\ . = we rewrite also effective 

. . W J mm 

potential Q5.22) in terms of variables y, A 



U, 



eff 



2D+4 



Ri 



.D-2 



(A 2 



di 

E 

i=0 



Pi (A 4+dl+1 + X dl - 1 ) 



Ay D + 
(5.34) 



The further analysis consists in a compatibility consid- 
eration of minimum conditions ( |5.30| ) and ( [5.32] ) with 
properties of the polynomial w(\), expressions like 
( |j.20 ) defining y D ~ 2 and y D = y D ~ 2 y 2 as functions 



crfi +3 
par 

i = 0, ...,di| and the constraint U e ff\ ■ =0. As 

" J C J J Imin 

result we will get a first crude division of Kp^ 3 in 
stability-domains allowing the existence of minima of 
the effective potential U e f f and forbidden regions cor- 
responding to unstable internal space configurations. 

After these general considerations we turn now to 
a more concrete analysis. 



of A on the parameter-space Rp^ 3 = {(Pi) A, Pi) 



5.2. Noncoinciding scale factors 

(A^l),i?i,A^0 

In this case we have to consider roots of type I ( |5. 19 ). 
We start from the polynomial u>(A) knowing that sta- 
ble internal space-configurations correspond to real 
projective coordinates < A < oo . So we have to test 
subpolynomials w±(\) on the existence of such roots. 
The high degree deg\[w±(\)] = 2(di+2)(di + l) > 24, 
di > 2 (because of nonvanishing curvature of the 
factor-spaces Mi , M2 ) allows only an analysis by tech- 
niques of number theory the theory of ideals 
of commutative rings | fl8l or, for general parameter- 
configurations, numerical tests. In the latter case the 
number of effective test-parameters can be reduced by 
introduction of new coordinates in parameter-space 



par 
par 




2Adi A dl+1 2P 



ftdl+2 ' 

(5.35) 



(for Pq ^ 0, po = 1 ; in the opposite case Po can be 
replaced by any nonzero Pi ) . Polynomials w± (A) = 
transform then to 

io±(A) = ±P? 1+2 ((di + l)P ) dl+1 u; ± (A) = 0, where 



w ± (X) := L d ^ +1 ± (-) dl+1 X &o 1+2 







L 3 :=L 3 {P = l;P 1 =p 1 ,...,P dl = Pdl ) (5.36) 

6 := b (P = l;P 1 =p 1 ,... ,P dl = p dl ). 
Test are easy to perform with programs like mathe- 

MATICA or MAPLE. 

As a second step we have to consider minimum con- 



ditions (5.30). Using (5.20) we substitute 
D-2 2 ^o . d _ ( d i + !)^3 



IJ 



y 



Adi 



(5.37) 



into (5.27) and transform (5.30) to the following equiv- 
alent inequalities 



2(di + l)(di + 2)Qi(A)+Q 2 (A) > 
> (d!+2)(A 2 + l)6 (A) 

[2(d 1 + 2)Q 1 (A)-(A 2 + l)fe (A)]x 
x[g 2 (A)-(A 2 + l)5 (A)] > 
>(d 1 + l)(A 2 -l) 2 5 2 (A) 

with notations 

QiW :=Eto^(A 4+dl+l + A dl -' ; ) 



(5.38) 



00(A) 



hW ■= E-io P S + 2 ) 2 (A 4+dl+i + A dl "0 



(5.39) 
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Stability-domains in parameter-space Kp,^ 3 , corre- 
sponding to minima of the effective pot ential are given 
as intersections of domains defined by ( 5.3S ) with do- 
mains which allow the existence of physical relevant 
roots of w±(X) — 0. So numerical tests on minima 
are easy to perform. If we additionally assume that 
U e ff\ . — then the class of possible st ability do- 
mains narrows considerably. Substitution of ( 5.37 ) into 
( 5.34 ) transforms this constraint to 

(A 2 + 1)6 (A) = (di+2)Qi(A) 

and inequalities ( |5.38| ) to 

di(A 2 + l)6 (A)+Q 2 (A) >0 

(A 2 + 1)[Q 2 (A)-(A 2 + 1)6 (A)]> 
>(di+l)(A a -l) 2 &o(A)>0. 



(5.40) 



(5.41) 



5.3. Coinciding scale factors (A = 1), R\, A 7^ 

In this case we have to consider roots of both types. 



Let us start with type I roots. Relations (5.37) read 
now 

Oa(l) 



.D-2 



y 



y D 

with 

Qi(i) 



Ri ' 
(di + l) 
Adi 



»=0 



[Q 2 (1)-2Q 1 (1)] 



(5.42) 



(5.43) 



Qa(l) = [(A 2 + 1)6 (A)] A=1 = 2]Tp i (* + 2) 2 

i=0 

so that at the extremum point the sign relation 
sign(J^ i= ]_o Pi) = sign(i?i) must hold and the corre- 
sponding scale factor is simply given as 



Vo 



rl \ 1 /(- D - 2 ) 



The minimum conditions (5.3q) reduce now to 



2(di+2)£P; >£P;(i + 2) 

i=0 i=0 



and for U e 

(rf 1 + 2)^P, = ^P,(i + 2) 2 >0. 



J ef f\ 
CJ J I min 

di 



— even to 



(5.44) 



(5.45) 



i=0 



i=0 



Substitution of (5.42) into (5.27) shows that for A = 1 
we have 



= y- 2D+2 (d 1 



^xxUef f\ x=1 — 9 xy U e ff\ x=1 



1) [2(di + 2)Qi(l)-Q 2 (l) 



(5.46) 



This implies det(A c ) — and according to ( 5.29 ) the 
eigenvalues W( c )i,2 of the Hessian A c ( |5.28 ) are given 



W( c )i = 2d yy U e ffi x=v w {c)2 
The condition uv c )i > is equivalent to the inequality 
(5.44). We note that the degeneracy ( pT46| ), det(A c ) = 
0, wt c )2 = is caused by the special and highly sym- 
metric form of the effective potential for the model with 
two identical factor-spaces. It is only in a secondary 
way related to the coinciding scale factors (x = y) 
for A = 1. In Appendix C we illustrate this fact with 
a simple counter-example, i.e., with a potential that 
gives W( c ji 2 > for a minimum at x = y. 

We turn now to the ro ots o f type II given by the 
polynomial h+{X = l,y) ( 5.21 ). From the structure 
of I 2 +(A = 1) immediately follows: 

1. Because I 2 +(A = 1) contains only terms with 
even degree in y , there exist no real roots — and hence 
no extrema of the effective potential U e t / — for pa- 
rameter combinations with: 



0. 



(5.47) 




sign(A) ^ sign(i?i 



(5.48) 



2. For arbitrary parameters A, i? 1; A := YliLo Pi 
roots of ^2+(A = 1) can be found by analytical meth- 
ods up to dimensions d\ < 2 performing a substitu- 
tion z := y 2 and using standard techniques for poly- 
nomials of degree deg 2 /2+(A = 1) < 4. Because of 
Ri ^ <^> di > 2 such considerations are restricted 
to the case di = 2 . 

3. There exist no general mathematical methods to 
obtain roots of polynomials with degree deg 2 /2+(A = 
1) > 4 and arbitrary coefficients analytically. For spe- 
cial restricted classes of coefficients techniques of num- 
ber theory Q are applicable. We do not use such tech- 
niques in the present paper. For polynomials ^2+(A = 
1 ) and dimensions dim M% — dim M 2 = d\ > 2 this 
implies that arbitrary parameter sets should be ana- 
lyzed numerically or parameters A, R±, A should 
be fine tuned — chosen ad hoc in such a way that 
h+{\ = 1) = is fulfilled. 

In the following we derive a necessary condition for 
the existence of a minimum of the effective potential 
with fine-tuned parameters. Using the ansatz 



Adi 
di + l 



(5.49) 



i=0 



equation ( |5.21| ) reduces to 



(a 1 -R 1 +4a 2 )y^- z = 0. (5.50) 

Without loss of generality we choose 02 as free param- 
eter, and hence a\ = R% — 4(72 , so that from relations 



(5.49) 



y^ 2 



A 

0-2 ' 



Vo 
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and (5.33) minimum condition (5.32) reads 



d yy U eff 



(5.52) 



4y - 2D+2 (d 1 + l) 



4(di + 2) - ^ 

^2 



A > 



or 



(2fl--)fp i >0, 
ff2 £S 



(5.53) 



We see that there exists a critical value er c = ^ which 
separates stability-do main s with different signs of A. 
From ~ 2 > and (|5.49|) follows sign(A) = sign(cr 2 ) 
so that 



A = Ef=o^<0 



CT 2 > cr c , cr 2 > 



(5.54) 



> (72 > cr c 



To complete our considerations of the degenerate case 
(A = 1), Ri, A 7^ w e deri ve the constraint ^e//| min = 
0. By use of (5.34) and (5.51) this is easily done to 
yield er 2 = -Ri /D — 2a c . So the constraint fixes the 
free parameter <y% . Remembering that according to 
our temporary notation y := a 2 = e 13 the value yo 
defines the scale factor of the internal spaces at the 
minimum position of the effective potential, we get now 



for the fine-tuning conditions (5.49) 



A 



A = 



QP-2)fli 

Da tc)2 

flf (£>-2) 



R\a 



D-2 
(c)2 



D 



D-2 



D D A D-2 



(5.55) 



— the well-known conditions widely used in literature 
[ p7[ . From ( 5.54 ), ( 5.55 ) and a( c ) 2 > we see that for 
<j 2 = Ri/D = 2a c the stability-domain in parameter- 
space R^i^" 3 is narrowed to the sector 



A 



par 
di 

E 

i=0 



Pi > 0, Ri > 0, A > 0. 



(5.56) 



5.4. Vanishing curvature-scalars 

(iZi = 0), A^0 



For vanishing curvature scalars equations (5.16) reduce 
to 



h+ = -2{D-2)a Q {X)-{D-4)ky 1 
h- = -2(A 2 -1)6 (A)=0. (6) 



(a) 
(5.57) 



Extrema of the effective potential are given by roots of 
I2- = with scale factors defined as 



,D 



2(di + l)oo(A) _ 2(di + l)Qi(A) 



Ad! 



(5.58) 



Substitution of ( 5.58 ) into minimum-conditions ( 5.3C ) 
yields the following inequalities 

Qi(A)>0, Q 2 (A)>0, Qi(A) + g 2 (A) >0 



8(di + l)(d 1 +2)Q 1 (A)Q 2 (A) > 
> (4di+5) 2 (A 2 -l) 2 ^(A). (5.59) 
From ( |5.58| ) and ( |5.59| ) we see that for even D positive 



y are only allowed when the bare cosmological constant 
A is negative: A < . 

As in the case of nonvanishing curvature scalars 
so also roots of 7 2 _ = split into two classes. For 
nondegenerate physical relevant configurations (A 7^ 1) 
the corresponding A^ must satisfy equation 

di i+l 

b {\) = ]T P t (2 + i)X dl -' 1 a2j = °- ( 5 - 60 ) 

i=0 j=0 

For A > this is only possible when there exist Pj 
with different signs. 

In the case of degenerate configurations (A = 1) 
equation J 2 _ = is trivially satisfied and the scale 
factor at the minimum of the effective potential given 
by 



„_ 4(di + l)Eto^ 

Vo — 



Adi 

with additional condition 



1] > 0. 



(5.61) 



(5.62) 



From inequalities (5.59) and (5.62) immediately follows 
that effective potentials with parameters (Po < 0, . . . , 
Pd 1 < 0) are not stable. 

5.5. Vanishing curvature scalars and 
vanishing cosmological constants 

(Ri =0, A = 0) 



In this case equations (5.16) contain only the projective 
coordinate A = y/x 

I 2+ = -2(D-2)oo(A) = (a) 

(5.63) 

h- = -2(A 2 -1)6 (A) = 0. (6) 

Corresponding physical configurations are possible for 
domains in parameter space given by 

di 

a = i, J2 Pi = ( 5 - 64 ) 



i=0 



A tU, i? A [a (A), 60(A)] =0. (5.65) 

Minima of the effective potential are localized at lines 
{Xi = y/x} and must be stabilized by additional terms. 
Otherwise we get an unstable "run-away" minimum of 
the potential. 
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5.6. Generalization to n— scale-factor models 

The analytical methods used in the above consider- 
ations on stability conditions of internal space con- 
figurations with two scale factors can be extended to 
configurations with 3 and more scale factors by tech- 
niques of the theory of commutative rin gs . In this 
case constraints, similar to polynomial (5.8) w(X), fol- 
low from resultant systems on homogeneous polyno- 
mials. We note that in master equations ( |5.5| ) we can 
pass from affine coordinates {x, y} to projective co- 
ordinates {X, Y, Z | x — X/Z 7 y = Y/Z} and trans- 
form polynomials I\± to homogeneous polynomials in 
{X, Y, Z} so that these generalizations are imme- 
diately to perform. A deeper insight in cxtrcmum 
conditions can be gained by means of algebraic ge- 
ometry jy|. Polynomials I\± define two algebraic 
curves on the {x, y}— plane and solutions of system 
( |5.5|) I\±{x,y) — correspond to intersection-points 
of these curves. For n scale factors extremum condi- 
tions {d ai U e ff = 0}™ =1 would result in n polynomials 
Ii(xi, . . . , x n ) = defining n algebraic varieties on 
W 1 . The sets of solutions of system {Ii(xi, ... ,x n ) = 
0}™ =1 , or equivalently, the intersection points of the 
corresponding algebraic varieties, define the extremum 
points of U e f f . 



6. Conclusion 

This paper is devoted to the problem of stable com- 
pactification of internal spaces in multidimensional cos- 
mological models. This is one of the most impor- 
tant problems in multidimensional cosmology because 
via stable compactification of the internal dimensions 
near Planck length we can explain unobservability of 
extra dimensions. With the help of dimensional re- 
duction we obtain an effective Do -dimensional (usu- 
ally Do = 4) theory in the Brans-Dicke and Einstein 
frames. The Einstein frame is considered here as the 
physical one. Stable compactification is achieved here 
due to the Casimir effect which is induced by the non- 
trivial topology of the space-time. The calculation of 
the Casimir effect in the case of more than one scale 
factors is a very complicated problem. That is the rea- 
son for proposing a Casimir-like ansatz for the energy 
density of the massless scalar field fluctuations. In this 
ansatz all internal factors are included on equal foot- 
ing. The corresponding equation has correct Casimir 
dimension: [cm] D and gives correct transition to the 
one-scale-factor limit. Stable configurations with re- 
spect to the internal scale factor excitations are found 
in the cases of one and two internal spaces. 
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Appendix A: Casimir effect 

It is well known that boundaries or nontrivial topol- 
ogy of the space-time induce a vacuum polarization 
of quantized fields due to changes in the spectrum of 
the vacuum fluctuations relative to a corresponding 
unbounded open flat-space model. This phenomenon 
is known as vacuum Casimir effect and in our case 
arises due to the compactness of the internal spaces 
Mi (i = 0, . . . , n) . In the general case of a quantum 
field living at finite temperature in a multidimensional 
cosmological model, the Casimir free energy due to the 
presence of compactified internal factor-spaces is de- 
fined as 



F c = F[{a i }^ 1 ]-F[{a i ^cx>}t i ] 



(A.l) 



where F [{<!,}?_]] and F [{at — > oo}" =1 ] are the free 
energy for a compactified (finite-scale-factor) model 
and a corresponding decompactified (flat) model, re- 
spectively. It is clear that this definition is only 
meaningful, when supplemented with a regularization 
method yielding finite results for the Casimir energy. 

In what follows we briefly discuss some regulariza- 
tion techniques for the Casimir energy of a massive 
scalar field in a universe with metric (|2.l|) and give 
some arguments in favour of an improved version of 
the Casimir-like potential, proposed in Ej|, |p4| . 

Let us consider a massive (with the mass M ) scalar 
field $ with arbit rary coupling to gravity in a universe 



with the metric (2.1). The Klein-Gordon- Fock equa- 



tion for the conformally transformed scalar field 
$ = $a n (1 - do)/2 



(A.2) 



reads [E8 



+ [M 2 al + £do(do - l)k + 



(A.3) 



n dj(dj—l)ki 



$ = o, 



where the overdot denotes the derivative with re- 
spect to the conformal time r\ (we choose the time 
gauge: e 7 = e 13 = do), £ is the coupling constant, 
Mi (i = 0, . . . , n) are the spaces of constant curvature: 
R [ g (i)] = kidi{di - 1), h = ±1,0, and A [gW] are 
the Laplacc-Bcltrami operators on Mj : 



A 



1 







VW\ dxr ' 



g ^) mn __ ]{AA) 



To get the pure Casimir effect without d ynam ic effects 
admixture we supposed in the equation ( A.3 ) all scale 
factors to be static, what implies that the space-time 
itself becomes ultrastatic. 
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We denote by x 1 the collective spatial coordinates 
of the i-th space and make for the scalar field the prod- 
uct ansatz 



^ = g{v)Y(x°)...Y(x n ) 



(A.5) 



Further, we assume that the scalar field on the cosmo- 
logical background is situated in the eigenstates of the 
Laplace-Beltrami operators 



A 



(i) 



Yi = -nfY, , i = 0, 



(A.6) 



so that (A.3) can be rewritten as 

-£di(di 



9 



M 2 



E 

1=0 



Thus the physical frequency squared of the scalar field 
reads 



M 2 



t=0 



+ £,di(di - l)ki 



(A.7) 



The free energy of the scalar field at temperature 
k B T = 1/0 is now given as & pi 



F 



= ^ In [1 -exp(-0tkJ fi 



.)] 



E< 



(A. 



where J is a collective index for all quantum num- 
bers. The eigenvalues — nf depending on the topology 
of the spaces M, and the boundary conditions imposed 
on Yi can be expressed in terms of the energy level 
quantum numbers rij with degeneracy Pi(rii). For ex- 



ample, n 
S 1 ■ ^ 

?2. ^2 



(" 



;/- = n(n + 1) (n 



0,±1,±2,...) , 

0,1,2,. 



p = 1 for 
) , p = 2n + 1 

for S 2 ; n 2 = n 2 - 1 (n = 1, 2, 3, . . .) , p = n 2 for 
S 3 and n 2 =n{n + 3) , p= (n + l)(n + 2)(2n + 3)/6 
for 5 4 . Our further consideration we restrict, for sim- 
plicity, to the case of the zero-temperatur e (v acuum) 
Casimir effect. Setting T = in equation (A..8) we get 
for the vacuum free energy 

2 E • • ' E P °( n °) • ■ -Pn{nn)uJn a ,...,n n = 

tlq n n 



F = 



J 



(A.9) 



so that the Casimir free energy according to (A.l) can 
be rewritten as 



dJu>„ 



(A.10) 



The internal energy density of the fluctuations is found 
from the relation 



1 d(0F c 
V 80 



1 



-,F r 



(A.ll) 



where V = n"=o ^ = n"=o a ? ^ s ^ ne s P ace volume 
of the universe and 



IH= f d di y\J\9^\, i = 0,...,n. 



(A.12) 



Pressure of fluctuations and Casimir energy- momentum 
tensor read, respectively, 



Pi = - 



1 



dF n 



1 



dF n 



Vdi da; 



i = 0, 



,n . 



and 



T 



M 
N 



diag 



-p, Pq, ■ ■ ■ , Po, 



p 



(A.13) 



(A.14) 



V 



do 



Let us now turn to the consideration of regularization 
techniques that are necessary for a meaningful defini- 
tion of the Casimir energy jAT| ), flATCl ). Most ap- 
propriate for the regularization of multiple sums like 
( A.9 ) proved zeta- function technique at the one hand 
(for an extended review see (ItJ and cites therein), 
and the multiple use of the Abcl-Plana summation 
formula on the other hand jl6). For simplicity, we 
perform our explicit consideration for a manifold M 
consisting of Ricci-flat (toroidal) factor-spaces M = 
K x S 1 x S 1 x • • • x S 1 . In this case the physical scalar 
field frequency reads 



/ 



LU 2 j = M 2 



in 2 



\ 



'01 



L 0d 



a 2 



\ external space internal 



(A.15) 



and cummulative index and degeneracy are given as 
J = {"oi, ■ ■ • ,n Mo ,m, ■ ■ ■ ,n n }, n x, ...,n n e Z and 
p.j{J) := Poi(«oi) ■ ■ -Pn{n n ) = 1. Because of p.j{J) = 
1 the free energy ( [A.9] ) can be expressed in terms of 
analytically continued inhomogeneous Epstein zeta- 
functions 

Z^(s; b N ) := 

:= (M 2 + b 1 7i\ + --- + b N n 2 N )~ s (A.16) 

ni : - ■ ■ ,un 

and 



E™{s;b x ,...,b N ) 



J2 (M 2 + bml + ■ ■ ■ + b N n 2 N ) s (A.17) 



directly fl7| 



n oi , ■ ■ ■ : n n = — oo 
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jV-1 

2 ^ 

J=0 



A/ 



E E N-j 
{ki,...,k s } 



S ^ 2' 



(A.18) 



where {k\ 1 . . . ,kj} denotes a combination of j indices 
- whose corresponding terms are absent — out of 
N = do+n possible ones and we have included the pref- 



actor 47r into the scale factors ai/2-ir 



for conve- 



nience. (The absent terms a^ 2 result from the zeros 
of the summation indices = when we change the 
summation sets from Uqi, . . . , n n € Z to Uqi, . . . , n n G 
{1,2,3, . . .} passing from Zff to Eff .) Multiple iter- 
ative use of the exact formula [B3| 



E%(s;b u ...,b N ) 



2\j h T(s) 



V 



r(s)" 1 ^fc!(167r) s 



fe/2 k 

L6^n [( 



2s - 1 



...,b N ) 

M + 

) 2 - (2j 



E 



rii,...,njv=l 

2 , u „2 , | t 2 \ -(«+*) /2 



x (Af 2 + b 2 n\ H h b N n 2 N ) 

2tt 



x exp 



(M 2 + 6 2 riJ 



(A.19) 



with start condition for the iteration 



<( S ;1) 



Af- 



/?Fr(a- 



r( s ) 



-M 



-2s+l 



2tt s M- s+1 / 2 



E 



pS-^K^/^TrpM) (A.20) 



would allow to express the Casimir energy in terms 
of exponential functions and modified Bessel functions 
K v of the second kind. With these formulae at hand 
an analysis of the asymptotic behavior of the Casimir 
energy can be performed. For the case of two scale 
factors with a 2 < a\ the corresponding Epstein zeta 
functions after an additional regularization have been 
estimated in Ref. [B4| to yield 



(A.21) 




, :exp(-2^) [1 + 0(1CT 3 )] 

2V«l a 2 V "2/ 



and 



ni ,n2 = l 





M 2 



(A.22) 



M 7T , 1 [M Mai \ 

= a 1 a 2 M s + — J x 

4 6 y An V a 2 Aira 2 ) 

x exp (-2nMa 2 ) [l + 0(1(T 3 )] . 

To our knowledge, for three and more scale factors 
similar estimates are performed up to now with less 
precission (see e.g. Jl^|). The calculations become 
more complicated, when one tries to develop analyti- 
cal approaches for compactified non-Ricci-flat internal 
factor-spaces with more then one scale factor. In this 
case the nontrivial degeneracy factor pj(J) ^ 1 in the 
sums of type (A.9) prevents a direct use of Epstein 
zeta function methods. Examples for an analytical 
circumvention of this problem by use of resummation 
and Mellin transformation are given in |l7j for the case 
of different one-scale-factor-spaces. Generalizations to 
multi-scale-factor models will need additional efforts in 
future. 

For completeness we make now some brief remarks 
on the regularization method based on the Abel-Plana 
summation formula IpOl, Bill, [£32] 



£/(») = ^ + / /<"W" 

71=0 



' / — r — -dv (A.23) 

exp (Ims) — 1 



for a function f(n) satisfying some conditions in the 
complex plane C . This method was develope d mai nly 
in Refs. fill , ||l6| and consists according to ( A.lOj ) in 

oo 

a subtraction of divergent terms containing J f(n)dn. 

o 

It yields the Casimir energy as finite expression. We 
applied this method, e.g., to get low-temperature cor- 
rections to the vacuum Casimir effect caused by scalar 
fields living, for example, in a closed Friedmann uni- 
verse p2[ and on a manifold with three-dimensional 
torus topology M = R x S 1 x S 1 x S 1 || (for a 
consideration of the high temperature regime with the 



help of zeta funct ion techniques see 
plication of ( A.23 ) to the triple sum ( A.K 
for a massless scalar field [E9| 



o 



Direct ap- 
gives, e.g., 



OO 

+ J dl 

o 



drll t!->0 
b a 



<h'[[ -;-,0 

c a j 
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dvl\ t;-,0 ) + / dvl(-;-,0 



-4 / / dn dm 

/ m 2^„i 

oo _ _ oo 



(v m n 
dvl\-; —, - 

\a o c 



-2 / dn 





E 

Z = /-5 5 



i/ I n 
b ' a' c 

\ 



E 



r v m n 
la' 6 ' c 



/ 



-4E 



Z,m=0 



dW 



i/ Z m 



CP. ^ ,(A.24) 



where 



1) \(A.25) 



a, 6, c are the scale factors corresponding to the wrapped 
factor-spaces of the three-torus, and CP. means terms 
obtained by cyclic permutations {a, I; b, m; c, n} , for 
exa mple, 6— — > c- — » a^ . Similar to the derivation 



of ( A. 21 ) the symmetri c ent ering of the factor-space 
scale factors a,b,c into ( A. 24 ) will be recovered explic- 
itly only at a final stage of an asymptotic expansion. 
For some considerations on a possible use of an Abel- 
Plana summation formula generalized to C" we refer 
to Appendix B. 

From the fact that the scale factors for factor-spaces 
of the same topological type (the same dimension and 
curvature) up to higher order corrections enter into the 
physical frequency according to (A. 7) in the same way 
we conclude that there must exist a corresponding in- 
terchange symmetry in the expression for the Casimir 
energy. Asymptotical expansion (A. 21) for the toroidal 
two-scale-factor model shows this symmetry explicitly. 
Furthermore, from the one-scale-factor model that ef- 
fectively emerges in the case of coinciding scale factors 
a\ = ci2 = . . . = a n = a we know that the Casimir 
energy density in this case has the form 



P = Ce- D ? = ^ 



C_ 

7d 



(A.26) 



where D = l + do + X)" = i ^ 1S the dimension of the to- 
tal space-time and C is a dimensionless constant which 
depends strongly on the topology of the model. It is 
clear that the Casimir energy density has the same 
structure for models with only one scale factor jul, 




U (or all scale factors identically equal 
cto = gi = . . . = a n = a ) . Explicit 



calculations, formula ( A. 21 ) from |34| and asymptotic 
expansion 



90a 4 



C(3) 

2nab 3 6abc 2 



(A.27) 



derived in papers jig] , p6[ for the case c > b > a 
show that the Casimir energy density as function of 
the scale factors has an effective degree deg(p) = — D 
and in a serial expansion of p terms containing a^ D , 
D > D are compensated by factors containing other 
scale factors. A sufficiently simple function with this 
feature and which under decompactification of some of 
the factor-spaces yields the correct new degree can be 
constructed as 



•£<=o 



1 - ' E e fcoifci...fc„i E E • 

feo,...,fc„=0 f o =0£i=0 



E 



£„-l=0 



(A.28) 



where A 



(fco 



are dimensionless constants which 



depend on the topology of the model and eik...m = 
(— l) P £ife...m ■ Here £ik...m is the totally antisym- 
metric symbol (eoi...m = +1), P is the number of 
the permutations of the 01 ... n resulting in ik . . .m. 
\k1k2 ■ ■ ■ k n \ means summation is taken over k\ < &2 < 
•■• < k n . Form ally, the Eq. (A.28) follows from 



expression ( A. 21 ) if we generalize it to the case of 
M = R x T d ° x • • • x P rf ", where T dl are df-tori, and 
omit the exponential terms. 

Let Mo denote the external space and its scale fac- 
tor Oq — > 00 . Then equation (A.28) yields 



E £ fc 1 |fe 2 ...fc„i E ■ 

fci,...,fe„=l 5i=0 



(A.29) 



where A 



(fci) 



,, , = Af 1 ) , and D Q = d + 1. For 
example, in the case of the two internal spaces and 



a ^> ai,a2 Eq. (A.28) reads 



P = 



E 



A 



a 



(i) 

* af 0+ « 



di 



E^ 



(2) 



.(A.30) 



It is clear, that for identical spaces Mi = M 2 we have: 



di = d 2 , P [ff^ 1 )] = P [<?( 2 )] and 4 = A 



I (2) 
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Equations ( |A.28|) , (|A.29| ) and (|A.30[) have correct 



Casimir dimension: cm and give correct transition 



to the case of one scale factor ( A.26). This can be 
easily seen, for example, from Eq. ( A. 30 ) when one of 
the scale factors goes to infinity. In one's turn equation 
( A.3C ) may be obtained from equation ( A. 29 ) if n — 2 
internal scale factors go to infinity. Thus, in contrast 
to the expression (1.3) for the Casimir energy density 
proposed in the papers j2^| , (24|] in our case the energy 
density does not equal to zero if at least one of the 
scale factors is finite. 

In the case of a large external space: ao S> a\ , . . . , a n 
from the equations (A. 11), ( A.13| ) and (A.2E) we can 
easily derive the equation of state in the external space: 



Pa = ~P , 



(A.31) 



because here F c = Vp = (Uto Mi) e 7 V (P\ ■ ■ ■ , P») 
and 70 = Yh=o ^ > and 



2J diPi = D p. 



(A.32) 



Appendix B: Abel-Plana summation 
formula in C n 

Here we give some arguments for a generalization of 
the Abel-Plana summation formula ( A.23| ) to the case 
of several variables. For this purpose we first recall 



that (A. 23) is the result of a contour integration in C. 
Given a contour 7 = dA C C as boundary of a rect- 
angle A — {z e C I \Im(z)\ < R;a< Re{z) < b} one 
considers the Cauchy contour integral 



f{z)ctg(irz)dz. 



(B.l) 



If necessary the contour 7 should be deformed ap- 
propriately to pass around possible branch cuts of 
f(z). Adding § dA+ f(z)dz-§ dA _ f(z)dz, where A± 
< Im{z) < R 
> Im(z) > R 
of the rectangle A = A + U A_ located at the upper 
and the lower complex half-plane, we arrive at 



zeC 



; a < Re(z) < b > are the halfs 



E 

s£L A 



f(s) 



f(z)dz + 
f(z)dz 



7+ 



f(z)dz 

^ g— I'Kiz 



- 1 



(B.2) 



z e 7 



and L A = (a, 6)nZ. To get Eq. 



Here the integration contours are defined as j± 
< Im(z) 
> Im (z) 

( |A.23|) the conditions \f(iR + t)\ < ip(t)e alR ^,a < 2tt; 
f a <p(t)dt < M < 00 and ip(t) — > should be satis- 
fied. This means that for an application of the Abel- 
Plana summation formula to the calculation of the 



Casimir energy the physical frequency should be mul- 
tiplied by an appropriate cutoff function. At the end 
of the calculations this cutoff can be removed. 

Relation ( |B.2| ) is well suited for a generalization 



to a multi-contour integral in C™. Following |36|, |37 
we start with a function f(z) polyholomorphic in the 
domain S := Ai x A2 x • • • x A n C C" built from the 
rectangles Ai c C}. Then for f(z) holds the integral 
representation 



(B.3) 



(2m) r 



fit) 



as 



(in %7i) 



dt\ A • • • A dt r , 



where dS is defined as dS := dA\ x dA 2 x • • • x dA n . 
Making use of J2k=~oo T^k ~ 7TC ^s{ 7ra ) we can extend 
this representation to 



{fci,...,fe„}eis 



(B.4) 



If " 

— * f(z 1 , z n ) TT [7rctg(7rz i )] dz\ A • • • A dz n 

n ) JdS f = i 



(2ni 



with Ls — S R Z n — the n— dimensional lattice seg- 
ment contained in the domain S . Proceeding as in the 
case of one complex variable, a relation generalizing 
(B.2) can be obtained and convergence criteria can be 
established. The main task in considerations for func- 
tions f(z±, . . . , z n ) with branch cuts, as in the case of 
the Casimir energy, consists in nontrivial contour de- 
formations and in keeping control over signs and sign 
changes of the integrals when moving along the ori- 
ented contours. 

As in the calculations of the Casimir energy in Refs. 
fig] , fl6| ] step by step calculations with the help of (B.4) 
lead to expressions containing terms with some func- 
tions g of the form g(^-zi) . These functions lead 
to an aparent asymmetry between the scale factors et^ 
at intermedium stages of the calculation, as explicitly 
shown, e.g., in (A. 24) and |l7j] also. Nevertheless, final 
asymptotic expansions can be transformed to expres- 
sions symmetric in the scale factors. 



Appendix C: Potential with minimum 
of type if( c )i,2 > for coinciding scale 
factors 

In this Appendix we show with the help of a sim- 
ple example that a minimum of an effective poten- 
tial U e ff(x,y) at a point with coinciding scale factors 
(x = y) implies not necessarily a degenerate Hessian 
A c d^j) det(l c ) = 0,w (c)1 > 0,w (c)2 = . As po- 
tential we choose 



U eff (x, y) = f(x)g(y) + h(xy) 



(C.l) 
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with 

f(x) — a^x^ 1 + a+x 
9{y) = b-y- 1 + b+y 
h(xy) = c^(xy)^ 1 + c + xy. 



(C.2) 



This ansatz is not related to any of the physical poten- 
tials considered as approximation to the exact Casimir 
potential. But we will see that for some coefficient sets 
{a±,b±,c±} it provides a minimum with w^ 12 > 
for coinciding scale factors x = y > . Proceeding 
as in subsection 5.1. we start from the extremum con- 
dition d x U e ff — 0, dyUeff = and consider the 
equation system xd x U e ff ± yd y ll e ff = 0. For x = y 
this yields 



-(a_6_ + c_)y- J + (a + b + + c+)y^ = 

(-o_6+ + a+b^y- 1 = 
so that for y > we have 

, CL-b- + c_ 



y 



(C.3) 



(C.4) 



a + b + + c + 

and the second derivatives of the potential read 

dxxUeff = dyyU e ff = 2(a + b + + c+) + 2a_6 + y -2 

2(a+6+ +c+) - 2a_6+y" 2 . (C.5) 



d xyUeff 



Hence the eigenvalues of the Hessian A c are given as 

W (c)l,2 = dyyUeff ± dxyU e ff Or W( c)1 = 4(fl + 6 + + C+) 

and Wt c )i t 2 — ^a-b + y~ 2 . Choosing appropriate non- 
vanishing coefficients a± , b± , c± the minimum condi- 
tions if(c)i,2 > are easy to satisfy so that in the 
general case the potential U e f / can have a nondegen- 
erate minimum at a point with coinciding scale factors 
x = y. 
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